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Abstract 

Quantum fc-SAT is the problem of deciding whether there is a n-qubit state which is perpendicular 
to a set of vectors, each of which lies in the Hilbert space of fc qubits. Equivalently, the problem is to 
decide whether a particular type of local Hamiltonian has a ground state with zero energy. We consider 
random quantum fc-SAT formulas with n variables and m = an clauses, and ask at what value of a 
these formulas cease to be satisfiable. We show that the threshold for random quantum 3-SAT is at most 
3.594. For comparison, convincing arguments from statistical physics suggest that the classical 3-SAT 
threshold is a c « 4.267. For larger fc, we show that the quantum threshold is a constant factor smaller 
than the classical one. Our bounds work by determining the generic rank of the satisfying subspace for 
certain gadgets, and then using the technique of differential equations to analyze various algorithms that 
partition the hypergraph into a collection of these gadgets. 

1 Introduction 

In quantum fc-SAT [T], each clause corresponds to a projection operator on the Hilbert space of n qubits, 

C= (lfe-|u)(u|)®l„-fc. 

Here \v) is a vector in the 2 fe -dimensional Hilbert space of some fc-tuple of qubits, 1*. is the identity 
on that Hilbert space, and l n -k is the identity on the remaining qubits. A formula is a set of clauses 
<f> — {C'i, . . . , C m }- We say that <f> is satisfiable if there is a state \rp) which is perpendicular to all the 
forbidden vectors \v): in other words, if 

(i/)\Ci\ip) = 1 for all i. 

We call the subspace of such states \tjj) the satisfying subspace Kat- 

In addition to being the quantum analogue of a canonical NP-complete problem, quantum fc-SAT is 
an illustrative case of a fc-local Hamiltonian. In that case, V SB x is the subspace spanned by eigenstates 
of a Hamiltonian H = ^,AI — Ci) with zero energy. It was shown in [T| that the decision problem of 
whether a particular quantum fc-SAT formula is satisfiable is in P for k = 2, and is QMAj-complete for 
k > 4, where QMAj^ is the subclass of QMA where the probability of acceptance for a yes-instance is 1. 

We are also interested in the problem of determining the rank J? sat = dim V sa t of the satisfying 
subspace, or equivalently the degeneracy of the zero-energy ground states. Determining 7? sat is a nat- 
ural quantum analogue of a classical counting problem, namely determining the number of satisfying 
assignments of a fc-SAT formula. Classically, even for k — 2 this problem is #P-complete under Turing 
reductions [2J. In the quantum case, it is not obvious that finding _R sat is even in =ffP, since the satisfying 
states may be arbitrarily entangled and may have no succinct description. Indeed, it seems to us that 
one can define a natural quantum version of #P as the class of problems consisting of finding the rank 
of an eigenspace of a fc-local Hamiltonian, although we do not pursue this further here. 
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In the classical setting, a lively collaboration between computer scientists, statistical physicists, and 
mathematicians has grown up around the behavior of random fc-SAT formulas. These are constructed 
in the following way. In order to construct a random formula <j>(n, m) with n variables and m clauses, 
we first construct a random fc-uniform hypergraph with n vertices and m edges, by choosing m times, 
uniformly and with replacement, from the (?) possible fc-tuples of vertices. Then, for each edge, we 
choose uniformly from the 2 fc possible combinations of signs for those k literals. 

We are particularly interested in the sparse case, where m — an for some constant a. There is a 
conjectured phase transition, where these formulas go from satisfiable to unsatisfiable when a exceeds a 
critical threshold: 

Conjecture 1. For each k > 3, there is a constant a c such that 

I 1 if d < ~C. c%c 

lim Pr[<^(n, an) is satisfiable] = < c 

n^aa I if a > Q c . 

In the absence of a proof of this conjecture, one can prove statements of the form that <j> is unsatisfiable 
with high probability if a > a*, or satisfiable with high probability if a < a}. Then, assuming that a 
phase transition exists, a* and o? are upper and lower bounds on the threshold a c . The state of the art 
for classical 3-SAT is [1 H [5] 

3.52 <a c < 4.490 , 

although overwhelmingly convincing arguments from physics [B] indicate that 

a c « 4.267. 

In the quantum case, we can similarly define a random quantum fc-SAT formula 4>{n, m) as a random 
hypergraph, where for each edge we choose the forbidden vector \v) uniformly from the vectors of norm 
1 in the Hilbert space Cf k of those fc qubits. We can then conjecture an analogous phase transition at 
a critical density a*. Laumann et al. [7] showed that 

0.818... < a\ < Q , 

where 0.818... is the density at which the hypergraph has a nonempty 2-core with high probability. In 
this paper, we show that 

a% < 3.594 . 

Note that this upper bound is well below the accepted value of the classical 3-SAT threshold. We also 
show that for all k > 4, 

a q c < 2 k b, 

where b w 0.573. Since the classical threshold grows as 2 fe ln2 [8], this shows that the ratio a%/a c is 
strictly less than 1. 

In order to prove these results, we exploit the observation of Laumann et al. that, once the hypergraph 
G is fixed, -Rsat takes a generic value R^t with probability 1. One way to see this is to note that 
with probability 1, the components of the clause vectors are algebraically independent transcendentals. 
Then any subdeterminant of the matrix of forbidden vectors is zero if and only if it is zero when these 
components are replaced by indeterminates. Moreover, for any particular choice of the clause vectors \ v) 
we have _R sat > Rf^ti since this choice can only result in linear dependences among the forbidden vectors 
and thus increase the rank of the satisfying subspace. 

Our bounds work by partitioning random hypergraphs into gadgets for which we can compute 
exactly. In order to show that certain partitions exist, we use the technique of differential equations to 
analyze simple greedy algorithms. To our knowledge, this is the first time that differential equations 
have been used to prove upper bounds on satisfiability thresholds. 



2 The case k = 2 

As a warm-up, in this section we reproduce results of Laumann et al. [7] on quantum 2-SAT, determining 
iif™ for all multigraphs, and in particular determining for which multigraphs the corresponding formula 
is generically satisfiable. 
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Theorem 1. Let G be a connected multigraph with n vertices and m edges. If we form a quantum 
2- SAT formula by replacing each edge (i,j) with a clause forbidding a random vector \vij) G C® 2 , then 
its generic rank R 1 , 



sat 




(G is a tree), 

,,m, ■ * uj ui — il (G is a cycle or a tree with a double edge), , . 

= > 1 ™-a (G consists of a triple edge), ( ' ^ 

Proof. Let V = {1, . . . ,n} be the vertices of G and E be its edges. Let <f> = {|iiij)}(i,j)e.E be a fixed 
instance of 2-SAT defined on G. Here \vij) £ C2 <8>C2 is a forbidden state associated with the edge 

Let O be an invertible local operator, or ILO — that is, O = ® iev Oi where all the O; are invertible, 
but not necessarily unitary. Define a new instance 

O • (j) = {Oi® Oj |«ij)}(i,j)ea;. 

We claim that <j) and O ■ (f> have the same rank. Indeed, a state \tjj) is a satisfying assignment for iff a 
state (O^) -1 \i[>) is a satisfying assignment for O ■ <f>. 

If <j> is a generic instance, all states are entangled. Let T be any spanning tree of G. Then there 
exists an ILO O that maps all forbidden states on the edges of T to singlets r^(|0, 1) — ] 1 , 0) ) [1]. The 
operator O maps forbidden states on edges (i,j) ^ T to some new forbidden states which are still generic 
(although their new distribution might not be uniform). Thus it suffices to compute Rf™ f° r instances 
<j> such that all edges of T are singlets and all other edges are generic states. 

Let (fitroe be the restriction of <j> onto the tree T (all clauses £ T are removed). Any satisfying 

assignment of </>t fe e is invariant under transpositions of any pair of qubits (i,j) £ T, and thus invariant 
under any permutation of qubits. Obviously, the converse is also true. Thus the satisfying subspace of 
</Vcc is exactly the totally symmetric subspace S n C C® n , which has dimension dimS n = n + 1. 

Suppose m = n — 1. Then G is a tree, <j> = <j>tree, and thus iif™ = n+1. 

Suppose m > n. Let \vij) be any clause of </>\<?!>trcc. Since the satisfying assignments of cf> span some 
subspace of 5 n , the choice of i and j doesn't matter — applying the clause \vij) to any pair of qubits gives 
the same rank. Let us apply all clauses of <j>\<t>trcc to the pair of qubits 1, 2. There are m — n + 2 forbidden 
states on this pair of qubits: the singlet from 0tree, and m — (n — 1) forbidden states from 0\(/> tre e. 

If m > n + 1 then there are at least 3 forbidden states on qubits 1, 2. These completely fix a state of 
qubits 1,2, say |wi,2). In the generic case |wi,2) is entangled. If n > 3, the monogamy of entanglement 
implies that \u)i t 2) cannot be symmetrically extended to n qubits. In that case, there are no satisfying 
assignments and Rf^t — 0- If n = 2 then |a>i,2) is the unique satisfying assignment, and -Rfa" = !• 

It remains to consider the case m = n, where G contains a single cycle or is a tree with a double 
edge. Now qubits 1, 2 have two forbidden states: the singlet and some state | t/^i , 2 ) - We can get an upper 
bound on by choosing IV'1,2) adversarially, for example, ^1,2) = |0, 1) + |1, 0). In this case we can 
look for satisfying assignments with a fixed number of Is, since all clauses commute with the particle 
number operator Y17=i F° r any number of particles < m < n there is only one symmetric state 

\S m ) — the uniform superposition of all binary strings with Hamming weight m. One can easily check 
that \Sm) is orthogonal to |0, 1) + |1, 0) iff m = or m = n. Thus there are two satisfying assignments: 
|0® n ) and |l® n ). This proves that 7? s g ^ n < 2. 

To show that Rf^t ^ 2, for any 1^1,2) we can try to construct a satisfying assignment |y>®") for some 
\tp) € C2. Without loss of generality IV'1,2} is symmetric, since otherwise it is some linear combination 
of the singlet and a symmetric state and we can redefine the clause. In the generic case 1*01,2) is also 
entangled. Consider the following proposition: 

Proposition 1. For any symmetric entangled state \if>) £ C2 ® C2 there exist two linearly independent 
states \ifi) , \ € C2 such that 

(ip\ip®ip) = (l>\<p'® <p') =0. (2.2) 

Proof. Let Aij — (%p\i,j) be the 2x2 complex matrix corresponding to \ifj). We are promised that A is 
symmetric, A T = A, and non-singular. Using Gaussian elimination, for symmetric matrices one can find 
an invertible complex matrix O such that 0^40 T = 1. This is equivalent to 

(O® 0)\i>) = 1 0,0) + |1,1> . (2.3) 
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Now we can choose 

|^)>') = (O t )- 1 (|0)± l |l)). (2.4) 

□ 

Thus = 2, and the proof is complete. □ 

Now suppose that we form a random multigraph with n vertices and m = an edges by choosing 
uniformly with replacement from the (™) possible edges. It is well known that if a < 1/2, then with high 
probability every connected component has at most one cycle, or one double edge, but never both — while 
if a > 1/2, then with high probability there is a giant connected component with multiple cycles. Thus, as 
already shown in [7], random quantum 2-SAT has a phase transition from satisfiability to unsatisfiability 
at a = 1/2. 

On the other hand, the classical 2-SAT transition is at a — 1 (see e.g. [TO]). This is a good illustration 
of the fact that the generic quantum problem is much more constrained. 



3 Expected gadget projectors 

The next lemma generalizes a result of Laumann et al. [7J , which showed that addiing a fc-clause reduces 
iif™ by a factor of 1 — 2 _fe . Our argument is somewhat simpler. 

Lemma 2. Let G and H be hypergraphs with n and t < n vertices respectively. Let G U H denote the 
hypergraph resulting from, adding a copy of H, on some subset ofG's vertices, to G. Then 

hsz(g u h) < 2- t R^(H)R^(G) . 

Proof. Let Hh be the projection operator onto the satisfying subspace of H , viewed as a subspace of 
Cf n . Now consider its expectation Ells, taken over the choice of clause vectors \v). Since each \v) is 
chosen uniformly from the sphere in Cf k , and since the uniform measure is invariant under any rotation 
of a single qubit, E YIh commutes with any one-qubit unitary operator affecting a vertex in H . Since YIh 
acts as the identity on the other n — t vertices, it commutes with one-qubit unitary operators on them 
as well. 

with any one-qubit operator affecting a vertex in H . Since YIh acts as the identity on the other n — t 
vertices, it commutes with one-qubit operators on them as well. 

Thus II_ff commutes with any Pauli operator. Since these form a basis for the full matrix algebra 
acting on Cf n , it follows that Erin must be a scalar. Since 

rkriH = r-'R^iH) 

holds with probability 1 (where the factor of 2 n_ * comes from being able to set the other qubits of G 
arbitrarily), it also holds in expectation. Thus 

trEITH = EtrR H = 2" -t i^"(#) , 

and therefore 

En H = 2- t Ri::(H)i. 

Now suppose the clause vectors of G are in general position. For any choice of clause vectors on H 
we have 

i?LT(G ufl)< rkn GuH < tr(n G n H n G ) . 

(This follows from the fact that ABA is positive whenever A and B axe projection operators.) This is 
also true in expectation over the clause vectors of H , so 

RSftGUH) < Etr(n G n H n G ) = 2- i 7?f a c t n (//)trn G = 2- t R!::(H)Ri::(G) , 

completing the proof. □ 
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Figure 1: The (6, 3)-sunflower. 



As pointed out in [7], if we take If to be a single clause for which = 2 fc — 1, this shows that for 
a random formula with n variables and m = an clauses we have 

fl£T < 2"(1 - 2~ fe ) m = [2(1 - 2" fe ) Q ] " . 

If a > log g(/7 2 « 5.191 then this bound is exponentially small, showing that such formulas are unsatis- 
fiable and placing an upper bound on the critical threshold. A similar argument applies in the classical 
case. However, in the next sections we will show that in the quantum case, we can prove much stronger 
upper bounds by computing R^t f° r larger gadgets. 



4 Two handy gadgets 

In this section we compute the generic rank exactly for two families of hypertrees. We will use these 
calculations to derive our upper bounds on the critical threshold. 



4.1 The sunflower 

Consider the (d, k) -sunflower, the fc-uniform hypergraph consisting of n clauses (edges), each pair sharing 
a common "center" vertex z. Specifically, the graph is defined over the 1 + d(k — 1) vertices 

{z} U {x\ I 1 < j < d, 1 < i < k - 1} , 

and contains the d clauses 

Cj = {z} U {x\ I 1 < i < k - 1} . 

See Fig. [1] for an example. 

Lemma 3. Let S(d,k) denote the generic rank of the (d,k)- sunflower. Then 

S(d,k)=2(2 k - 1 -l) d ( ¥ ^ + l 
Proof. Decompose the Hilbert space of the (d, fc)-sunflower as 

where 

Ho = C 2 , Hi = --- = H d = (C 2 f k ~ 1 . 



5 



Figure 2: The (1, 2, 3) -nosegay. 



Here Ho describes the central qubit and Jij describes the other k — 1 qubits on the jth petal for 1 < j < d. 
Let Cj denote the clause on the jth petal, and let \vj) G TLo ®TLj be its forbidden state. Clearly \vj) has 
at most two non-zero Schmidt coefficients, so we can always choose a unitary operator Uj acting on TLj 
such that 

\v j ) = (I ®U j ) \uj) ® |0® fe " 2 ) . (4.1) 

Here \uj) is some entangled state between the central qubit and the first qubit of the jth petal. Let us 
refer to the remaining k — 2 qubits of the jth petal that are projected onto |0) as ancillas. Clearly Uj 
does not change the rank, so without loss of generality we can assume that Uj =1 for all j. 

If we ignore the ancillas, then we get an instance of quantum 2-SAT on a star graph with d edges. 
By Theorem [1] its generic rank is d + 2. Clearly ancillas of the jth petal can be ignored iff they are all 
set to the state |0). Let us say that such a petal is active. Otherwise, if at least one ancilla of the jth 
petal is |1), it becomes inactive because the corresponding clause is already satisfied. An inactive petal 
contributes a factor of 2 k ~ 1 — 2 to the rank. For a fixed subset A of active petals, the generic rank is 

R(A) = (a + 2){2 k - 1 -2) d - a , a=\A\. (4.2) 

Therefore, the generic rank for the (d, fc)-sunflower is 

R= J2 R(A) = ^( d \a + 2)(2 k - 1 -2) d - a =2(2 k - 1 -l) d ( ¥ ^- + l) . (4.3) 

AC{l,...,d} a=0 V / ^ ' 

This completes the proof. □ 

4.2 The nosegay 

For another example, consider a nosegay, as shown in Fig. [2] It consists of a single edge, where each of 
its vertices has some number of additional edges attached to it. 

3-uniform nosegays Let us momentarily restrict our attention to the case k = 3. Then an (a, b, c)- 
nosegay has a, b, and c additional edges, for a total ofa + 6 + c+ l edges and 3 + 2(a + b + c) vertices. 

Lemma 4. Let R( a ,b,c) denote the generic rank of the 3-uniform (a,b,c)-nosegay. Then 
R(a,b,o) = 3 a+b+c ~ 3 [(a + 6)(b + 6) (c + 6) - (a + 3) (b + 3) (c + 3)] 



G 



Proof. Let us label the qubits of the central triangle by 1, 2 and 3. These qubits have a, b and c hanging 
triangles attached to them respectively. Each triangle represents a generic forbidden 3-qubit state. 

Let us also define the [a, b, c]-nosegay: it coincides with the (a, b, c)-nosegay except that each hanging 
triangle is replaced by a hanging edge. Accordingly, the [a, b, c]-nosegay has n — 3 + a + b + c qubits. 
Each hanging edge represents a generic forbidden 2-qubit state while the central triangle represents a 
generic forbidden 3-qubit state. 

Let R( ai b,c) and R[ a ,b,c] be the generic ranks of the (a, b, c)-nosegay and the [a, b, c]-nosegay. Repeating 
the arguments used to compute the genetic rank of sunflowers in Lemma [3J we get 

In the rest of the section we prove that 

R[a,b,c] = (a + 2)(b + 2)(c + 2) - (o + 1)(6 + l)(c + 1) (4.5) 

which after simple algebra yields 

R( a>b ,c) = 3 a+6+c - 3 [(a + 6)(6 + 6)(c + 6) - (a + 3)(6 + 3)(c + 3)] . (4.6) 

We shall start from using the symmetry of the [a, b, c]-nosegay to bring the forbidden states into a 
canonical form such that the forbidden state associated with the central triangle is 

Ka,s> = -^(|0,0,0)-|l,l,l» (4.7) 

while the forbidden state associated with any hanging edge is the singlet, 

K-> = -L(|0,l)-|l,0». (4.8) 



We claim that any set of generic forbidden states can be mapped to the ones defined in (|4.7[1 . (|4~i 
by applying invertible local operators (ILO) to every qubit. Indeed, it was shown by Diir, Vidal, and 
Cirac [11] that a generic 3-qubit state is ILO-equivalent to the GHZ state which is in turn ILO-equivalent 
to the state |ui,2,3) defined in (14. 7[) . Note that applying an ILO at this step maps the forbidden states on 
the hanging edges to some entangled 2-qubit states. We can convert them to singlets by applying proper 
ILO to the free end of every hanging edge. As we argued in Theorem [T] the dimension of the satisfying 
subspace is invariant under ILO, so it suffices to compute R[ a ,b,c] with forbidden states defined by (|4.7[1 
and gU). 

Our arguments will rely on the fact that the canonical forbidden states define an instance of sto- 
quastic 3-SAT studied in [T2]. An instance of stoquastic 3-SAT is defined by a family of projectors 
<j> = (IIi, • • • ,II m ) which have real non-negative matrix elements in the computational basis, such that 
every projector acts on at most 3 qubits. A state is a satisfying assignment for <j> iff II a — for 
every a = 1, . . . , m. In our case the instance (j> is defined by a family of projectors 

ni,2,3 = 1- |«1,2,3)(«1,2,3|, n 4 ,j = 1 - \Vij)(Vij\ (4.9) 

where 111,2,3 acts on the central triangle and EEjj acts on every hanging edge A direct inspection 

shows that matrix elements of the above projectors have the following properties: 

• Any off-diagonal matrix element belongs to the set {0, 1/2} 

• Any diagonal matrix element belongs to the set {1, 1/2} 

A family of projectors with such properties defines an instance cf> of simplified stoquastic 3-SAT [121 
Section 6.3]. In particular, it was shown in [T2] that the number of satisfying assignments of <f> is equal 
to the number of connected components of a graph G — (V,E), where V = {0, l} n and (x,y) £ E iff 
there exists a projector II a £ (j> such that (x|II a |j/) = 1/2. (A satisfying assignment associated with a 
connected component V a C V is the uniform superposition of all vertices in V a , see [T2] for details.) It 
remains to count connected components in the graph G associated with the [a, b, c]-nosegay. 



7 



Let us partition the qubits of the nosegay into 3 disjoint subsets A, B, C such that qubit 1 and the free 
ends of all hanging edges attached to it form A, qubit 2 and the free ends of all hanging edges attached 
to it form B, and the remaining qubits form C. By definition, 

\A\=n a = a + l, \B\=n b =b + l, \C\=n c = c+l. (4.10) 

Using the projectors that live on the hanging edges we conclude that (x, y) G E whenever the restrictions 
of x and y onto any of the subsets A, B, C have the same Hamming weight. Leaving out the projector 
111,2,3 temporarily we can thus label the connected components of G by triples of integers 

(a,/3,7), 0<«<n„, < P < n b , < 7 < n c , (4.11) 

where for any vertex x £ V we define a, /3 and 7 as the Hamming weight of x\a, x\b, and x\c respectively. 
Adding the projector Hi,2,3 adds extra edges to the graph G, which glue together some components 
according to the following rules: 

(q j( 8, 7 )~ (a + l,/3+ 1,7 + 1), (a, 0, 7) ~ (a - 1, - 1, 7 - 1). (4.12) 

Thus the number of connected components of G is the same as the number of diagonals parallel to the 
axis (1, 1, 1) in a cube of size [0, n a ] x [0, rib] x [0, n c ]. This yields (|4.5|l and completes the proof. □ 



/c-uniform nosegays A fc-uniform nosegay is determined by a vector d = (di, . . . , dk) of nonnegative 
integers. It is the k- uniform hypergraph given by a single edge A — {cti , . . . , aj,}, the ith vertex of which 
is incident upon di other "hanging" edges B}, B^ with Bf = {a;} U {$(£) | 1 < I < k - 1}. These 
"hanging" edges intersect the central edge at a unique vertex, and are otherwise all disjoint. We refer 
to such a graph as a d-nosegay. In the next lemma we determine the rank of its satisfying space when 
adorned with separable clause vectors, which is an upper bound on its rank in the generic case. 

Lemma 5. Let N(d) denote the rank of the satisfying subspace of the d-nosegay when adorned with 
separable clause vectors in general position. Then 



N(d) = H(2 fc 



n(*+2(2 fc - i -i))-n(^+( 2fc_ 



Proof. As the clause vectors are separable, we may introduce a basis for the Hilbert spaces (copies of C 2 ) 
associated with the vertices fi{ (£) so that the clause vector associated with B\ has the form \a J i ) Cg> lO}* -1 . 
(Here the first factor in this tensor product is associated with the vertex aj.) As with the sunflower, we 
may expand a satisfying vector \v) according to this basis: 

\v) = Y, l b > ® K> > 

b 

where |b) = Cg> I&2) 55 • • • is a basis vector in the tensor product of the Hilbert spaces associated with 
the vertices b\ (£) and Vb lies in the Hilbert space Ha associated with the center edge. 

We write b h B\ when b t — 1 for one of the indices associated with B\ . Should b \f B\ , observe that 
we may expand «b in a Schmidt decomposition \ Vs ) ® where the \v 3 ) lie in the Hilbert space 

associated with 014, and, considering that the \w a ) are orthogonal and that \v) satisfies the clause B\, 
conclude that each (v a ,a{) = 0. It follows that \vb) has the form Io0 ® \w) (where Ia0 is orthogonal 
to I ) ) . As the |a^) are in general position, then, for each i we must have \{B? \ h h B?}\ < 1. 
Additionally, observe that if, for each i, we have b h Bf for some j then Vb is completely determined 
(and, in fact, separable), and cannot satisfy the clause A. In general, writing k — £ = \{Bj | b h B\ }|, 
we find that orthogonality with the vector associated with A precisely constrains Vb to a subspace of 
TLa of dimension 2 —1. Evidently, this expresses the satisfying subspace as an orthogonal direct sum of 
subspaces of total dimension 

£ (n^ 2 " 1 - i ) d '~ i ) ■ (n ( 2fc_i - if ) ■ ( 2fcHl1 - ■ 

IC{1,...,A} J J 

equal to the expression in the statement of the lemma. □ 
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5 Upper bound on the critical density 



In this section we present two upper bounds on the critical threshold. The first one is weaker but simpler. 

Theorem 6. Let H be a random ^-uniform hypergraph with n vertices and m = an edges. If a > 3.894, 
then with high probability the corresponding quantum 3- SAT problem is unsatisfiable. 

Proof. Let If be a fc-uniform hypergraph with n vertices and m edges, corresponding to a formula with 
m clauses. If we can partition H into a set of sunflowers, where there are nd sunflowers of each degree 
d, then by Lemmas [2] and [3] the generic rank of its satisfying subspace is bounded by 

^ a c t n <2™n^ s{d ' k) x 



d=l 



I 2 l+<*(fc-l) 



*S((>-^)' 



d 



2 k - 2 



+1)1 (5.1) 



2n n 



d=l 



3\ d (d 



47 V6 +1 



where in the last line we set k — 3. 

Clearly (|5.1|l is minimized if we have a small number of sunflowers of high degree. Ideally, we would 
like to characterize the best possible such partition. For now, we content ourselves with the partition 
resulting from the following simple algorithm: at each step, choose a random vertex, declare it and its 
edges to be a sunflower, and remove them from the graph. 

We can carry out this algorithm in continuous time, by assigning each vertex a uniformly random 
index t £ [0, 1] and removing vertices in the order of decreasing t. In that case, by the time we remove 
a vertex v with degree t, its sunflower includes those edges whose other vertices all have index less than 
t. Since this is true of each of its clauses independently with probability and since the original 

degree distribution of v is Poisson with mean ka, the degree of v's sunflower at the moment when it is 
removed is Poisson-distributed with mean fccrf fe_1 . Integrating over t, the expected number of vertices 
whose sunflowers have degree d is nad, where 



, . at " 1 (kat k - 1 ) d , 
a d = I -dt 



f 

Jo 



(5.2) 



1 e~ 3at (3at 2 ) d J 
d! 

_ r(rf+l/2) -r(rf+l/2,3a) 
2 % /3^d! 

where T(a, z) = J°° x a ~ 1 e~ x dx is the incomplete Gamma function. 

We then upper bound iif™ by cutting off the product above d max = 100, ignoring the effect of the 
tiny fraction of sunflowers of greater degree. Standard Azuma-type inequalities tell us that, with high 
probability, the number of sunflowers of degree d is adn+o(n) for all d < rf max . Also, with high probability 
there are less than log n pairs of edges which share more than one vertex. Their neighborhoods consist of 
sunflowers with two petals stuck together. We claim that such a sunflower has lower rank than a normal 
one, but in any case pretending that these petals are not stuck together only changes the rank by a 
constant factor, and the effect of logn such steps changes the rank by poly(n). So, with high probability, 



«S < p»iy(») 



2 ni(I) (f+> 



and therefore 



lim -InRf" < ln2+ Y~ a d ( dln-+lnf - + 1 
n^oo n / — ' V 4 \ 6 



If we set a = 3.894, we find that this limit is —1.372 x 10 , so is exponentially small. Thus with high 
probability in H , -Rf™ = with probability 1 in the clause vectors, and the formula is unsatisfiable. □ 
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Next, we improve this result by partitioning the graph into nosegays instead of sunflowers. Although 
the analysis is slightly harder, the algorithm is equally simple. 

Theorem 7. Let H be a random 3-uniform hypergraph with n vertices and m = an edges. If a > 3.594, 
then with high probability the corresponding quantum 3- SAT problem is unsatisfiable. 

Proof. At each step we choose a uniformly random edge, declare it and the edges it shares a vertex with 
to be a nosegay, and remove them and its vertices from the hypergraph. The remaining hypergraph has 
3 fewer vertices. Moreover, if we condition on how many edges it has, it is uniformly random in the 
model where edges are chosen with replacement. This allows us to model this process with differential 
equations [13] . 

If t is the number of steps we have taken so far, then there are n — 3t remaining vertices. Let m 
denote the number of remaining edges. Its expected change on each step is 

„. . , H 9m 
E[Am] = -1 - 



3t 

Now we write m = /m and t = rn, and rescale this to give a differential equation: 



^ = -1-^-. (5.3) 
dr 1 — 3r 



Changing variables to the fraction v = 1 — 3r of vertices remaining, this is 



^4 + ^. (5.4) 
du 3 v 



With the initial condition — a, the solution to this is 



6 

This becomes zero when vo vertices are left, where 



^{v) = ^-({6a + l)v 2 -1) . (5.5) 



vq = 



v 7 ^ 



At that point, there are no edges left, and the algorithm stops. 

With high probability, for any v > the number of edges remaining when there are vn vertices left 
is m(y) = \i(ii)n + o(n). Similar to the proof of Theorem HJ] summing over the (1 — vo)n/3 + o(n) steps 
of the algorithm then gives 



lim - In Ri:: <ln2+i 

n— >oo n J 



hi 



R 



(a,6,c) 



23 + 2(a+6+c) 



dv 



where R^ a ,b,c) is given by Lemma [4] and where a, b, and c are chosen according to independent Poisson 
distributions with mean 3/i/V. 

If we set q = 3.594 and upper bound the expectation over a, b, and c by ignoring terms where 
any of them is greater than 50, then evaluating the resulting integral numerically we find that limit 
is —1.601 x 10 -4 . Again i s exponentially small, so these formulas are unsatisfiable with high 

probability. □ 

There are a number of potential ways to improve this result. First, we can achieve a better partition 
of the graph into gadgets by prioritizing high-degree vertices. Analogous to Q3], we can analyze the 
resulting partition using a system of coupled differential equations, using the configuration model to 
keep track of the random graph conditioned on its degree distribution. For the sunflower, this gives a 
bound of 3.689 — a significant improvement over Theorem|SJ but not as good as Theorem]?] We have not 
attempted a partition into nosegays that prioritizes high-degree clauses. 

Second,we have no obligation to consider partitions into sunflowers or nosegays that can be found 
in polynomial time. We could also use non-algorithmic proofs that a desirable partition exists. These 
algorithms simply happen to be both efficient and easy to analyze. 

Thirdly, we could use notions of local maximality which have been successful in the classical case, 
but it is not obvious how to apply these in the quantum setting. When is an entangled satisfying state 
locally maximal? 
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6 A upper bound for general k 



In this section we use our sunflowers to prove an upper bound on the quantum fc-SAT threshold for 
general fc. We have made no attempt to optimize this bound beyond the simplest possible argument, but 
it establishes that the quantum threshold is strictly less than the classical one for all k > 6. 

Theorem 8. Let b m 0.573 be the unique positive root of the equation In 2 — 2b + \n(b + 1) =0. Then for 
all k > 3, if a > 2 b then with high probability the corresponding quantum k-SAT problem is unsatisfiable. 

Proof. First we rewrite (|5.ip as follows: 

«st< 2- (i- ^r) m n (^2+1 

since ndd = m. Treating the product as a harmonic mean over the vertices, bounding it as an 
arithmetic mean, and using the fact that the mean degree of a sunflower is E[d] = a then gives 



- In R»S < In 2 + a In I 1 - ) - L 



' ii '2^2 +1 



<m2 + aln(l-^ T j+ln(^ + lj . (6.1) 
Rearranging, applying the Taylor series of ln(l — x), and setting a = 2 k b, we have 
-laR^ <ln2+(a-l)ln(l-2 1 " fc ) + In (V fc a + 1 
< In 2 - (a - l)(2^ fe + 2 1 ~ 2fc ) + In (V fc a + 1 
= ln2-26 + ln(6+l) + 2 1 " fc 



2~ k a + 1 



= 2 1 ~ fe (1-6 — + 2" fe 

V 6+1 

< for all k > 3 , 

since 1 - b- 1/(6+ 1) < -1/8. □ 

In contrast, the classical fc-SAT threshold is known to be (l + o(l))2 fe ln2 [8]. Since b < ln2, it follows 
that the quantum threshold is less than the classical one for sufficiently large k. In fact, explicit lower 
bounds on the classical threshold for finite k from 8. are greater than the upper bounds on the quantum 
threshold obtained by setting (|6.1|) to zero for k > 6. Using the approach of Theorem [6] improves this to 
k > 5. 



7 Open questions 

We close with several open questions. 

• What is the computational complexity of determining the generic rank of a hypergraph? It would 
be surprising if it were not at least #P-hard, but it is not obvious that it is in #P. On the other 
hand, we are not aware of any proof that it is even NP-hard. 

• Can we prove a lower bound on the satisfiability threshold which is greater than the density at 
which a random graph contains a non-vanishing 2-core? In particular, is there a phase where 
random formulas are satisfiable, but all satisfying states are entangled? 

• Assuming that the quantum fc-SAT threshold exists, is al proportional to 2 k ? If so, what is 
b = \im k a q c /2 k ? The value of b given in Theorem [5] is almost certainly an overestimate. Note 
that our lower bounds based on the existence of the 2-core as given in [7 actually decrease as k 
increases — for instance, for k = 4, 5, and 6 we have the lower bounds 0.772, 0.701, and 0.637. So, 
at present, we do not even know that a.% grows without bound. 
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